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ABSTRACT 


The  number  of  containers  handled  by  container  terminals  has  increased 
significantly  over  the  last  fifty  years  and  has  stimulated  researchers  to  improve 
storage  yard  operations.  Container  handling  and  crane  deployment  are  two  major 
yard  operations  that  can  impact  the  performance  of  a  whole  container  terminal. 
This  thesis  establishes  an  Integer  Linear  Program  (ILP)  to  combine  container 
handling  and  yard  crane  deployment  for  Rubber  Tired  Gantry  Cranes  (RTG). 
Using  real  world  data,  we  test  the  ILP  for  two  different  yard  sizes.  We  find  the 
resulting  ILPs  difficult  to  solve  directly.  In  order  to  decrease  the  computation  time, 
we  apply  a  cascade  method  that  solves  the  problem  as  a  sequence  of  restricted 
subproblems.  Each  subproblem  is  restricted  to  a  sequence  of  containers  and  the 
output  of  each  subproblem  provides  an  input  to  the  next  subproblem.  This 
method  provides  better  solutions  than  the  solution  that  we  get  by  solving  the 
problem  directly.  The  cascade  method  also  decreases  the  computation  time 
significantly.  The  results  demonstrate  the  ability  to  combine  container  handling 
and  yard  crane  deployment  in  a  single  model  and  they  verify  that  the  cascade 
method  works  well  with  the  ILP. 
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EXECUTIVE  SUMMARY 


The  rapid  growth  of  containerization  significantly  increased  the  number  of 
containers  and  the  competition  among  seaport  container  terminals.  As  a  result  of 
competition  among  container  terminals,  managers  and  operators  began  to 
search  for  new  methods  to  improve  terminal  performance. 

Berth  time  of  a  vessel  is  usually  accepted  as  the  primary  measure  of 
performance  for  container  terminals,  and  this  time  depends  on  a  sequence  of 
operations  that  are  conducted  in  the  container  yard  where  containers  are  stored 
temporarily.  Loading  and  unloading  containers  and  assigning  yard  cranes  to 
specific  areas  in  the  yard  are  two  major  yard  operations  that  can  impact  the 
performance  of  a  whole  container  terminal,  unless  they  are  well  planned  and 
coordinated.  Although  there  have  been  various  studies  to  improve  these 
operations  separately,  this  thesis  is  the  first  to  establish  a  combined  model  for 
both  operations. 

In  this  thesis,  we  minimize  the  time  of  placing  a  container  which  is 
unloaded  from  a  vessel  (discharge)  in  the  container  yard  and  the  time  of  loading 
a  container  to  a  vessel  in  such  a  way  that  the  same  type  of  discharge  containers 
are  stored  together  and  the  total  workload  is  divided  fairly,  as  much  as  possible, 
among  Rubber  Tired  Gantry  Cranes  (RTG).  We  let  the  workload  of  an  RTG 
exceed  a  fair  limit  when  it  is  worth  paying  a  deviation  penalty.  We  define  the 
planning  horizon  in  terms  of  the  sequence  of  containers. 

We  formulate  the  problem  as  an  Integer  Linear  Program  (ILP).  The 
objective  function  seeks  to  minimize  the  operation  time  which  includes  transfer 
times  of  RTGs  between  areas,  average  time  that  an  RTG  spends  with 
discharging  or  loading  a  container,  time  to  carry  containers  from  the  vessel  to  a 
storage  area  or  vice  versa,  and  the  reshuffling  time  that  a  container  causes  if  it  is 
stored  with  different  type  of  containers. 
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Using  real  world  data,  we  test  the  ILP  for  two  different  yard  sizes.  We  take 
60  containers  as  a  base  case  for  our  study.  The  ILP  contains  about  1,400,000 
equations,  507,000  continuous  and  487,000  binary  variables  and  takes  more 
than  twenty-four  hours  to  provide  a  near  optimal  solution.  In  order  to  decrease 
the  computation  time  and  to  satisfy  memory  requirement,  we  apply  a  cascade 
method  that  solves  the  problem  as  a  sequence  of  restricted  subproblems.  Each 
subproblem  is  restricted  to  a  sequence  of  containers  and  the  output  of  each 
subproblem  provides  an  input  to  the  next  subproblem.  This  method  provides 
better  solutions  than  the  solution  that  we  get  by  solving  the  problem  directly.  The 
cascade  method  also  decreases  the  computation  time  significantly.  The  results 
demonstrate  the  ability  to  combine  container  handling  and  yard  crane 
deployment  in  a  single  model  and  they  verify  that  cascade  method  works  well 
with  the  ILP. 


I.  INTRODUCTION 


A.  BACKGROUND 

The  number  of  containers  handled  by  container  terminals  has  increased 
significantly,  as  international  sea  freight  transportation  has  grown  rapidly  over  the 
last  fifty  years.  As  a  result  of  competition  among  container  terminals,  managers 
and  operators  began  to  search  for  new  methods  to  improve  terminal 
performance.  The  berth  time  of  a  vessel  is  usually  accepted  as  the  primary 
measure  of  performance  for  container  terminals  and  this  time  depends  on  a 
sequence  of  operations  that  are  conducted  in  the  container  yard  where 
containers  are  stored  temporarily.  Loading  and  unloading  containers  and 
assigning  yard  cranes  to  specific  areas  in  the  yard  are  two  major  yard  operations 
that  can  impact  the  performance  of  a  whole  container  terminal,  unless  they  are 
well  planned  and  coordinated.  Although  there  have  been  various  studies  to 
improve  container  handling  and  yard  crane  deployment  (separately),  this  thesis  is 
the  first  to  establish  a  combined  model  for  both  operations. 

This  thesis  establishes  an  Integer  Linear  Program  (ILP)  to  combine 
container  handling  and  yard  crane  deployment  strategies.  We  choose  to  model 
the  deployment  of  Rubber  Tired  Gantry  Cranes  (RTG),  among  different  types  of 
yard  cranes,  because  RTGs  are  used  extensively  in  container  terminals  around 
the  world.  Our  ILP  finds  the  best  place  in  the  yard  under  given  constraints  for 
containers  that  are  going  to  be  unloaded  from  a  vessel;  it  also  picks  an 
appropriate  container  from  several  possible  locations  to  be  loaded.  We  find  the 
area  covered  by  each  RTG  that  minimizes  the  time  to  load  and  unload  containers 
by  minimizing  the  interblock  movements  of  the  RTGs,  as  well  as  balancing  the 
workload  among  them. 
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1. 


Container  Terminals 


Container  terminals  are  the  facilities  where  cargo  containers  are 
transshipped  between  different  transport  vehicles  for  onward  transportation.  A 
container  terminal  can  be  partitioned  into  three  areas:  quayside  (or  berth)  where 
containers  are  loaded  or  discharged  onto/from  a  vessel,  container  yard,  and  the 
landside  (gatehouse)  which  is  the  landside  entrance  of  a  terminal  where  the 
container  inspection  and  paperwork  is  done,  and  container-carrying  trains  and 
trucks  come  into  or  go  out  of  the  terminal  (Figure  1 ). 


Figure  1 .  Container  Terminal  Layout  and  Flow  of  Operations 

In  container  yards,  containers  are  stacked  in  areas  called  blocks.  A  block 
usually  consists  of  approximately  twenty  containers  in  length,  six  to  eight 
containers  in  width  and  four  to  six  containers  in  height.  The  length  of  a  block  is 
called  bay,  the  width  and  height  are  named  row  and  tier,  respectively.  Therefore, 
the  location  of  a  container  in  a  block  is  defined  by  its  bay,  row  and  tier  number. 
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2. 


Terminal  Equipment 


A  container  or  a  cargo  container  is  a  large,  typically  metal,  box  used  for 
transporting  goods  (Figure  2).  There  are  several  ways  to  classify  containers, 
namely  by  dimension,  by  weight,  and  by  content.  Of  these  classifications, 
dimension  is  the  most  common  and  TEU  (twenty-foot  equivalent  unit)  is  the 
standardized  unit  of  measure.  Forty-foot  and  forty-five-foot  long  containers  are 
counted  as  two  TEUs.  Dimensions  of  a  short  container  (so  called  standard 
container)  are  20x8x8.5  or  20x8x9.5  (in  feet),  whereas  dimensions  of  a  forty  foot 
long  (Two  TEU)  container  are  40x8x8.5  or  40x8x9.5. 


Figure  2.  40-foot  container  (Left)  and  20-foot  container  (Right)  (From:  Huynh  and 

Walton  [2005]) 


Containers  are  transported  by  trucks  or  Automated  Guided  Vehicles 
(AGV)  (Figure  3)  at  terminals.  Trucks  that  are  used  only  in  terminal  operations 
are  named  Terminal  Tractors  (Figure  4)  or  Internal  Trucks  (IT),  whereas  the 
trucks  that  transport  the  containers  into  terminals  or  out  of  terminals  are  called 
External  Trucks  (ET).  In  some  container  terminals,  Automated  Guided  Vehicles 
carry  containers  between  quay  cranes  and  the  container  yard  with  their  positions 
controlled  via  wires  or  transponders. 
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Figure  3.  Automated  Guided  Vehicle  (From:  Gottwald  Port  Technology  [2007]) 


Figure  4.  Terminal  Tractor  (From:  Kalmar  Industries  [2007]) 

There  are  different  types  of  cranes  used  to  load  and  unload  containers. 
Cranes  that  load  and  unload  containers  from/to  vessels  are  called  Quay  Cranes 
(QC)  (Figure  5).  These  cranes  play  a  major  role  in  the  performance  of  a 
container  terminal.  Other  types  of  cranes  are  usually  used  to  stack  containers  in 
container  yards.  The  three  main  types  of  these  cranes  are  Automated  Stacking 
Crane  (ASC)  (Figure  6),  Rail  Mounted  Gantry  Crane  (RMG)  (Figure  7),  and 
Rubber  Tired  Gantry  Crane  (RTG)  (Figure  8).  Of  these  three  types  of  yard 
cranes,  ASC  is  the  only  one  that  can  operate  without  an  operator.  ASCs  move  on 
rails  and  they  don’t  have  the  ability  to  move  between  blocks.  Because  RTGs 
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operate  on  rubber  tires,  they  are  much  more  versatile  than  ASCs  and  RMGs. 
While  the  movements  of  ASCs  and  RMGs  are  restricted  to  a  block  or  to  the 
blocks  on  the  same  lane,  respectively,  RTGs  can  move  between  blocks  in  a  yard, 
even  between  the  ones  that  are  located  on  different  lanes.  The  only  issue 
concerning  the  movements  of  RTGs  is  that  more  time  is  required  to  change  lanes 
than  to  change  blocks  within  the  same  lane.  When  an  RTG  is  transferred  to  a 
block  on  a  different  lane,  it  moves  out  of  the  block  to  the  open  space  at  the  end 
of  the  current  block,  makes  a  90-degree  turn,  moves  in  parallel  to  the  block 
width,  lines  to  the  correct  lane  of  an  adjacent  block,  makes  a  90-degree  turn 
again,  then  enters  the  block.  These  90-degree  turns  take  extra  time  that  delays 
not  only  the  operation  of  the  RTG,  but  also  the  traffic  that  flows  on  the  road  when 
an  RTG  is  turning. 
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Figure  6.  Automated  Stacking  Crane  (From:  Kalmar  Industries  [2007]) 


Figure  7.  Rail  Mounted  Gantry  Crane  (From:  Doosanheavy  Industries  and 

Construction  [2007]) 
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Figure  8.  Rubber  Tired  Gantry  Crane  (From:  Kalmar  Industries  [2007]) 


There  are  also  other  types  of  vehicles  and  equipment  in  container 
terminals  such  as  straddle  carriers,  forklifts,  reachstackers  and  trains.  However, 
these  are  not  considered  in  this  thesis. 

Containers  that  are  unloaded  from  a  vessel  and  are  sent  to  a  yard  to  be 
stacked  are  called  discharges  or  import  containers,  and  the  ones  that  are  sent 
from  a  container  yard  to  a  quay  crane  to  be  loaded  onto  vessel  are  called  loads 
or  export  containers. 

3.  Terminal  Operations 

Basic  operations  on  the  quayside  are  discharging  and  loading  ships. 
When  containers  are  discharged  from  the  vessel  by  quay  cranes,  ITs  and  AGVs 
carry  those  containers  to  the  container  yard,  and  then  carry  the  load  containers 
from  the  yard  to  vessels. 

There  are  several  operations  conducted  in  container  yards.  Yard  cranes 
place  discharge  containers  in  blocks,  pickup  the  load  containers,  and  load  them 
onto  ITs  or  AGVs  when  a  vessel  is  berthed  at  a  port.  These  cranes  also  load 
import  containers  onto  ETs  and  unload  export  containers  from  ETs  when  there  is 
no  waiting  discharge  or  load  move  from/to  the  vessel.  They  also  relocate  some 
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of  the  containers  to  facilitate  future  loads  when  they  are  idle.  Because  ship 
turnaround  time  is  one  of  the  major  measures  of  performance  of  terminals, 
quayside  operations  have  priority  over  yard  and  landside  operations. 

On  the  landside  at  the  gatehouses,  paperwork  is  done  for  the  entrance 
and  exit  of  ETs.  There  might  be  container-handling  operations  at  the  train  station 
if  there  is  a  railway  access  to  a  container  terminal. 

B.  ORGANIZATION  OF  THE  STUDY 

This  study  includes  five  chapters.  Chapter  II  provides  an  overview  of  some 
previous  studies  that  have  been  conducted  to  improve  container  terminal 
operations.  Chapter  III  presents  the  problem  that  is  the  source  of  this  study.  In 
the  same  chapter,  the  characteristics  of  the  model  that  combines  the  container 
handling  and  crane  scheduling  operations  are  explained.  Chapter  IV  describes 
the  test  data  provided  by  Navis  Lie  and  and  presents  results.  Chapter  V  includes 
conclusions  and  suggestions  for  further  studies. 
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II.  LITERATURE  REVIEW 


A.  INTRODUCTION 

Due  to  the  numerous  advantages  that  containerization  brought  to 
international  trade,  it  has  become  one  of  the  most  essential  parts  of  international 
sea  freight  transportation  over  the  last  fifty  years.  The  rapid  growth  of 
containerization  significantly  increased  the  number  of  containers  and  the 
competition  among  seaport  container  terminals  [Steenken  et  al.  2004], 

Berth  time  of  a  vessel  is  usually  accepted  as  a  measure  of  performance 
for  container  terminals  and  this  time  depends  on  a  sequence  of  operations  that 
are  conducted  in  the  yard.  Container  handling  and  yard  crane  scheduling  are  two 
major  yard  operations  that  can  easily  become  a  bottleneck  for  a  container 
terminal  unless  they  are  well  planned  and  coordinated.  Although  there  have  been 
various  studies  to  improve  container  handling  and  crane  scheduling  separately, 
this  thesis  is  the  first  to  establish  a  model  that  combines  both  operations. 

B.  PREVIOUS  STUDIES 

As  Steenken  et  al.  [2004]  report,  containers  are  designed  to  make 
international  transportation  of  goods  easier  by  using  a  unit-load  concept. 
Cordeau  et  al.  [2005]  state  that  containerization  requires  less  product  packaging, 
it  also  reduces  damages  and  yields  higher  productivity  during  various  handling 
processes. 

Zhang  et  al.  [2006]  report  that  containerization  has  grown  at  an  annual 
rate  of  nine  percent  recently.  They  also  state  that  it  is  expected  that  by  2010, 
ninety  percent  of  all  international  sea  freight  will  be  containerized.  Due  to  this 
high  rate  of  growth  and  the  high  cost  of  terminal  structure  and  equipment 
investments,  terminal  operators  have  been  searching  for  more  efficient  ways  of 
handling  containers  without  enlarging  the  footprint  of  terminals. 
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Although  all  terminal  operations  are  strongly  interrelated,  due  to  the  large 
number  of  decisions  and  the  multi-objective  nature  of  the  problem,  the 
uncertainty  and  the  complexity  of  the  decisions  result  in  a  focus  on  one  or  two 
specific  operations. 

Some  recent  research  areas  include  berth  allocation  (e.g.,  Cordeau  et  al. 
[2005]),  quay  crane  scheduling  (e.g.,  Park  and  Kim  [2003],  Kim  and  Park  [2004]), 
ship  planning  (e.g.,  Steenken  et  al.  [2001]),  automated  guided  vehicles 
scheduling  (e.g.,  Kim  and  Bae  [2004],  Rashidi  and  Tsang  [2005]),  container 
handling  (e.g.,  Kozan  and  Preston  [1999],  Alvarez  [2006]),  crane  routing  (e.g., 
Kim  and  Kim  [1999a]),  straddle  carrier  routing  (Kim  and  Kim  [1999b]),  crane 
deployment  (e.g.,  Chung  et  al.  [2002]),  automated  stacking  cranes  scheduling 
(e.g.,  Zyngiridis  [2005]),  and  classification  of  container  terminal  operations  (e.g., 
Steenken  et  al.  [2004]). 

Because  getting  real-time  solutions  is  essential  in  solving  container 
terminal  problems,  computational  complexity  of  the  problem  is  a  real  concern  for 
researchers.  Hence,  heuristic  optimization  techniques  have  been  proposed  by 
several  researchers  (e.g.,  Linn  and  Zhang  [2002]). 

Daganzo  [1990]  studies  a  queuing  problem  that  arises  at  multipurpose 
port  terminals.  He  presents  queuing  models  to  predict  the  stochastic 
characteristics  of  the  traffic  overflow  from  a  multipurpose  terminal  to  the  rest  of  a 
port. 

Castilho  and  Daganzo  [1993]  focus  on  container  import  operations  at 
container  terminals.  They  present  methods  for  measuring  the  required  amount  of 
handling  effort  based  on  two  different  strategies,  one  of  which  tries  to  keep  all 
stacks  the  same  size  and  the  other  that  segregates  containers  according  to 
arrival  time. 

Kozan  and  Preston  [1999]  use  genetic  algorithms  to  reduce  the  container 
handling/transfer  times  and  ships’  turnaround  times  by  speeding  up  handling 
operations.  They  report  that  a  scheduled  storage  policy  where  containers  are 
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stored  closer  to  the  berth  is  better  than  a  random  storage  policy.  They  also  state 
that  the  storage  area  fullness  does  not  have  a  significant  effect  on  either  the 
random  or  the  scheduled  storage  policy. 

Kim  and  Kim  [1999]  formulate  the  optimal  routing  of  a  single  yard  crane  as 
a  mixed  integer  program  (MIP),  to  minimize  the  container  handling  time  of  a 
transfer  crane  that  includes  the  setup  time  at  each  yard-bay  and  travel  time 
between  yard-bays. 

Zhang  et  al.  [2002]  formulate  the  deployment  of  RTGs  among  blocks  as  a 
MIP,  and  solve  it  by  Lagrangean  relaxation  to  determine  the  routes  of  crane 
movements  in  a  way  such  that  the  total  delayed  work  is  minimized.  It  is  assumed 
that  the  locations  of  the  discharge  containers  are  known  prior  to  the  operation. 

Linn  et  al.  [2003]  formulate  the  yard  crane  deployment  problem  as  a  MIP 
to  investigate  dynamic  crane  deployment  in  the  container  yard  on  a  shift-to-shift 
basis  and  solve  it  with  a  least  cost  heuristic  method.  Again,  it  is  assumed  that  the 
locations  of  the  discharge  containers  are  known  prior  to  the  operation. 

Chung  et  al.  [2002]  also  study  the  problem  of  scheduling  the  yard  crane 
movements.  They  formulate  the  problem  as  a  MIP  and  solve  it  by  using 
Lagrangean  decomposition.  They  also  report  a  method  named  piecewise-linear 
approximation  that  is  efficient  for  large-size  problems  and  again,  it  is  assumed 
that  the  locations  of  the  discharge  containers  are  known  prior  to  the  operation. 

Zyngiridis  [2005]  focuses  on  the  problem  of  automated  stacking  crane 
scheduling  and  develops  Integer  Linear  Programs  to  prescribe  routes  for  one  and 
two  equally  sized  ASCs.  He  finds  that  one  ASC  working  alone  over  four  hours 
requires  up  to  70%  more  time  than  two  ASCs  working  together  to  accomplish  the 
same  required  container  movements. 

This  thesis  differs  from  the  previous  studies  by  combining  the  container 
handling  and  yard  crane  deployment  strategies  to  improve  the  efficiency  of  yard 
management  in  container  terminals.  This  thesis  does  not  aim  to  provide  real  time 
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solutions  in  a  few  seconds.  Our  goal  is  to  provide  a  near-optimal  solution  in  a 
reasonable  time  that  can  be  a  benchmark  for  heuristic  solutions. 
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III.  MODEL 


A.  INTRODUCTION 

As  mentioned  in  Chapter  II,  two  of  the  major  measures  of  performance  for 
container  terminals  are  vessel  berth  time  and  quay  crane  rate.  However,  all 
terminal  operations  are  interrelated  and  a  delay  in  one  of  those  operations  may 
cause  a  delay  in  quay  operations  and  decrease  the  performance  of  a  container 
terminal.  Container  handling  and  crane  scheduling  are  two  of  these  yard 
operations  that  can  easily  cause  a  performance  bottleneck  for  a  container 
terminal  unless  they  are  well  planned  and  coordinated. 

Although  all  terminal  operations  are  interrelated,  the  uncertainty  and 
complexity  of  the  large  number  of  decisions  and  the  multi-objective  nature  of  the 
problems  prevent  us  from  formulating  all  terminal  operations  in  a  single  model. 
Therefore,  we  focus  on  two  major  operations  in  a  container  terminal.  One  of 
those  operations  is  distributing  discharge  containers  in  a  container  yard  and  the 
other  one  is  deployment  of  RTGs  to  the  areas.  When  a  container  is  discharged 
from  a  vessel  by  a  quay  crane,  it  is  loaded  onto  an  internal  truck  or  on  an  AGV  to 
be  carried  to  the  yard  where  it  is  going  to  be  stored  by  the  yard  cranes.  It  usually 
takes  a  few  minutes  to  carry  a  container  from  the  berth  to  the  yard,  depending  on 
the  distance  between  them  and  the  speed  of  the  carrying  vehicle.  There  is  an 
interval  called  push  rate  between  discharge  container  arrivals  to  the  yard.  Push 
rate  can  be  adjusted  to  prevent  congestion  in  the  yard.  However,  it  is  supposed 
to  be  as  small  as  possible  to  minimize  the  discharging  time  of  the  vessel. 

In  this  thesis,  we  minimize  the  time  of  discharge  and  load  operations  in 

such  a  way  that  the  same  type  of  discharge  containers  are  stored  together  and 

the  total  workload  is  divided  fairly,  as  much  as  possible,  among  RTGs.  We  let  the 

workload  of  an  RTG  exceed  the  limit  when  it  is  worth  paying  a  deviation  penalty. 

We  define  the  planning  horizon  in  terms  of  the  sequence  of  containers.  Our  aim 

is  to  keep  this  horizon  as  long  as  possible  but  in  real  operations,  conditions  in  the 
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yard  change  and  RTG  assignments  to  areas  must  be  reevaluated  periodically. 
Additionally  runtime  of  the  model  and  the  memory  of  the  computers  we  use  limit 
the  number  of  containers  that  can  be  handled. 

Regarding  the  above-mentioned  characteristics  of  the  problems,  we  have 
the  following  assumptions: 

•  The  sequence  of  the  movement  type  (either  discharge  or  load)  is 
known. 

•  Only  RTGs  are  used  for  stacking  containers  in  the  yard  and  their 
initial  positions  are  known. 

•  All  RTGs  are  the  same  type  and  size.  Therefore,  they  cannot  cross 
over  each  other. 

•  When  an  RTG  needs  to  be  transferred  from  one  block  to  another, 
the  delay  it  causes  in  the  yard  traffic  is  included  in  RTG  transfer 
time. 

Container  handling  includes  placing  discharge  containers  in  the  container 
yard  and  picking  up  the  load  containers  from  appropriate  areas  (among  the 
alternatives).  We  want  the  same  type  of  containers  to  be  stored  in  the  same 
area,  so  we  penalize  the  discharge  moves  that  do  not  store  similar  containers 
together.  We  also  want  to  ensure  that  a  load  container  is  obtained  from  a  specific 
area  where  that  type  of  container  is  already  stored,  so  we  provide  the  locations  of 
those  containers  in  the  data.  We  also  ensure  that  all  containers  are  handled. 

Crane  deployment  covers  the  deployment  of  RTGs  in  specific  blocks 
where  they  operate,  until  all  containers  are  handled.  We  want  each  RTG  to  be 
matched  with  one  or  more  areas.  Thus,  each  RTG  is  matched  with  the  area 
where  it  is  initially  located.  When  an  RTG  needs  to  change  its  area,  a  cost  is 
defined  in  crane  minutes  and  it  is  assumed  that  this  cost  includes  the  cost  of  the 
delay  of  the  yard  traffic  that  the  RTG  causes.  When  an  RTG  handles  a  container, 
it  cannot  be  assigned  to  handle  the  next  few  containers  in  the  sequence  because 
it  can  delay  yard  operations  while  waiting  for  the  RTG  to  be  available.  To  prevent 
RTG  congestions  in  areas,  we  only  let  one  RTG  be  matched  with  an  area. 
However,  an  RTG  can  be  matched  with  multiple  areas.  We  calculate  the 
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workload  of  an  RTG  by  adding  up  the  time  that  an  RTG  spends  while  it  is 
changing  areas  (transfer  time)  and  the  average  time  it  takes  to  discharge  or  to 
load  a  container.  Therefore,  the  workload  is  represented  by  crane  minutes. 

B.  FORMULATION 

We  formulate  the  problem  as  an  Integer  Linear  Program  (ILP).  The 
objective  function  seeks  to  minimize  the  operation  time  which  includes  transfer 
times  of  RTGs  between  areas,  average  time,  that  an  RTG  spends  with 
discharging  or  loading  a  container,  time  to  carry  containers  from  vessel  to  the 
storage  area  or  vice  versa,  and  the  reshuffling  time  that  a  container  causes  if  it  is 
stored  with  different  type  of  containers. 

1.  Indices 


t,t'  eT 
a,  a' 


d  e  D 
I  e  L 

9,9' 


containers  in  the  time  sequence  of  their  movement 
te{  1,2,...,T} 

area  where  containers  are  stored  or  picked  up  for  loading 
ae{1,2,...,A} 

containers  to  be  discharged  cfcF 

containers  to  be  loaded  /  c  T 


RTG 


ge{1,2 . Gj 


2.  Sets 


t'  e  Nt  : 

(a,g)  e  InitiaIG : 
a  e  Loc, 


Set  of  containers  that  an  RTG  cannot  handle  if  it  handles 
a  container  at  time  t. 

Set  of  area  crane  pairs  where  RTG  g  is  initially  located. 
Set  of  areas  where  load  container  /  can  be  obtained. 
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3. 


Scalars 


TOTALG 

AREAS 

MOVESt 


Total  number  of  RTGs  in  use 

Upper  limit  for  the  number  of  areas  that  an  RTG  can  cover 
Upper  limit  for  the  number  of  movements  of  an  RTG 
between  areas  up  to  container  t 


4.  Parameters 


Travel  a 

PenDCat 
Transfer^ 
Initial a,g 
ADTa 
ALTa 


Time  to  carry  a  container  from  a  vessel  to  an  area  or  vice 
versa. 

Penalty  of  discharging  container  t  into  area  a. 

Time  to  transfer  an  RTG  from  area  a'  to  area  a. 

1  if  RTG  g  starts  in  area  a. 

Average  discharge  time  for  area  a. 

Average  loading  time  for  area  a. 


Deltat 

CmlWLgt 

PEN 


Shows  how  much  the  workload  of  an  RTG  can  exceed  the 
average  workload. 

Cumulative  workload  of  an  RT G  up  to  container  t. 

Penalty  for  exceeding  the  workload  limit. 


5.  Positive  Variables 

WLg  :  Workload  of  an  RTG. 

Devgt  :  Deviation  from  the  sum  of  Average  WL  and  Delta. 

6.  Binary  Variables 

DCagt  :  1  if  container  t  is  moved  by  RTG  g  to  area  a. 
LDggt  :  1  if  container  t  is  moved  by  RTG  g  from  area  a. 
Yag  1  if  RTG  g  covers  area  a. 
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Zagt  1  if  RTG  g  has  control  over  area  a  at  time  t. 

/W,  ,  :  1  if  RTG  g  moves  from  area  a'  to  area  a  at  time  t. 

a  ,a,g,r  ^ 

7.  Objective  Function 

(1)  Min  X(  Travel, +  PenDCat)-DC.tt 

a,g,t 

+  2  7 ~raVela-LDa,g,t 

a,g,t 

+  2  TmnSfera  ,a  ' M a’,a,g,t 

a',a,g,t 

+  Y.PENDevgt 

g,t 


8.  Constraints 


(2) 

IDC,,,=1 

a.g 

Vd 

(3) 

Z  Z  LD..,J 

aeLoC/  g 

=  1  V/ 

(4) 

Z  Z  Dc.,,r  +  s  1  -ZDC,,,  + 

a  t'sNt 

a 

(5) 

DCa,g,t  ~  Za,g,t 

\/a,g,t>  1  and  teD 

(6) 

LDa,g,t  ~  Za,g,t 

Va,g,t>  1  and  teL 

(7) 

n,9=i 

\/{a,g)  e  InitiaIG 

(8) 

Zn^i 

g 

Va 

(9) 

AREAS  Vg 

a 

(10) 

z„.„  *  Y,„ 

Va,g,f  >  1 

(11) 

z.w=1 

V(a',g)  e  InitiaIG 

(12) 

Zz«=i 

a 

vg,t 

(13) 

^a,g, f+1  ~  Za,g,t  + 

X  M  Va,g,t  <  T 

a',a,g,t 

(14) 

M,  *  Za',9, 

a,a,g,t 

Vg,a',a(a'  *a),t 

17 


(15) 

(16) 


07)  IXa!?i<1 

a 

<18>  Z  ^  MOVES, 

a',a,t'<t 


Vg,a',a(a'  *a),t>  1 
Vg,a',a(a'  *  a),f  >  1 

V(a',g)e  InitiaIG 
Vg.f 


(19) 

(20) 


M,,  =  £  Transfer,, a  +  Z(^OT.  *  +  *  LO.J 


Z^V 


I  ^ 


<  — -  +Deltat  +  Dev  t  Vg,t 

G 


9.  The  Objective  Function  and  Constraint  Description 

The  objective  function  seeks  to  minimize  the  time  and  the  penalty  of 
placing  all  containers,  as  well  as  minimizing  the  time  for  transferring  RTGs 
between  areas.  Constraint  sets  (2)  and  (3)  ensure  each  discharge  and  load 
container  are  handled.  Constraint  set  (4)  prevents  a  container  to  be  matched  with 
an  unavailable  RTG.  Constraint  sets  (5)  and  (6)  require  that  when  a  container  is 
discharged  to  area  a  or  it  is  loaded  from  area  a  by  RTG  g  at  time  t,  for  RTG  g  to 
have  control  over  that  area  at  time  t.  Constraint  set  (7)  states  that  each  RTG  has 
to  be  assigned  to  the  area  where  it  is  initially  located.  Constraint  set  (8)  limits  an 
area  to  be  matched  only  with  one  RTG.  Constraint  set  (9)  limits  the  number  of 
areas  that  an  RTG  can  cover.  Constraint  set  (10)  ensures  an  RTG  can  have 
control  over  an  area  if  it  is  assigned  to  that  area.  Constraint  set  (11)  ensures 
each  RTG  is  in  its  initial  location  at  the  start.  Constraint  set  (12)  ensures  each 
RTG  can  be  only  in  one  area  at  time  t.  The  following  four  constraint  sets  regulate 
the  transfers  of  RTGs  between  areas:  Constraint  set  (13)  states  if  an  RTG  has 
control  over  area  a  at  time  t+1,  it  has  to  have  control  over  the  same  area  at  time  t 
or  it  has  to  be  transferred  to  that  area  from  some  area  a'  at  time  t.  Constraint  set 
(14)  states  that  if  an  RTG  is  transferred  from  area  a'  to  area  a  at  time  t,  then  it 
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has  to  have  control  over  area  a'  at  time  t.  Constraint  set  (15)  states  if  an  RTG  is 
transferred  from  area  a'  to  area  a  at  time  t,  then  it  has  to  have  control  over  area 
a  at  time  t+1.  Constraint  set  (16)  states  if  an  RTG  is  transferred  from  area  a'  to 
area  a  at  time  t,  then  it  has  to  have  control  over  area  a'  at  time  t  and  it  has  to 
have  control  over  area  a  at  time  t+1  as  well.  Constraint  set  (17)  states  that  if  RTG 
g  moves  from  area  a'  at  the  start,  it  has  to  be  located  in  that  area  initially. 
Constraint  set  (18)  limits  the  interblock  movements  of  each  RTG.  Constraint  set 
(19)  defines  the  workload  of  each  RTG.  Finally,  constraint  set  (20)  limits  the 
workload  of  an  RTG  at  time  t  to  the  average  workload  and  a  specific  amount 
called  Delta  at  time  t  plus  the  deviation. 
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IV.  COMPUTATIONAL  STUDY 


A.  TEST  DATA 

Navis  LLc  provided  a  data  set  consisting  of  over  6,000  containers,  half  of 
which  are  discharge  containers  and  the  other  half  are  load  containers  [Ayik  and 
Golbasi  2007],  Navis  LLc  reports  that  the  data  set  shows  a  daily  workload  of  a 
typical  container  terminal. 

We  are  given  the  container  names  and  the  locations  of  the  load 
containers.  A  load  container  may  be  obtained  from  several  different  areas.  We 
know  the  distances  between  quay  cranes  and  areas.  We  assume  that  ITs  carry 
containers  at  a  constant  speed.  Time  to  carry  a  container  from  a  quay  crane  to 
an  area  or  from  an  area  to  a  quay  crane  is  found  by  dividing  the  distances  by  the 
IT-speed.  There  are  costs  associated  with  RTG  movements;  not  only  for  inter¬ 
area  movements  but  also  for  the  in-area  movements.  There  are  also  costs 
associated  with  the  areas  where  a  container  is  discharged.  Cost  of  storing  a 
container  with  a  different  type  of  container  is  higher  than  storing  it  with  containers 
of  its  type. 

B.  ANALYSIS 

We  break  the  data  set  into  parts  so  that  we  can  get  a  solution  for  each 
small  data  set  in  a  reasonable  time  without  exceeding  the  memory  limits  of  our 
computers.  In  this  study,  we  take  one  hour  as  a  reasonable  time  to  get  a  solution. 
We  use  GAMS  [GAMS  2007]  Version  2.0.33.5  and  Cplex  Version  10.0.1  to  solve 
instances  using  a  desktop  computer  with  3.73  GHz  CPU  and  3.00  GB  of  RAM. 

We  consider  a  portion  of  the  container  terminal  and  two  quay  cranes 
discharging  and  loading  a  vessel  that  berths  at  a  port.  We  consider  two  quay 
cranes  just  to  simplify  how  we  obtain  a  list  of  containers.  The  list  of  containers 
could  easily  be  a  representative  of  more  than  two  quay  cranes.  There  are  six 
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RTGs  operating  in  the  yard  that  consists  of  twelve  lanes,  each  of  which  has  three 
blocks.  The  lanes  are  parallel  to  the  berth  (Figure  9).  Initial  positions  of  the  RTGs 
are  randomly  chosen.  It  is  assumed  that  only  the  sequence  of  the  discharge  and 
load  moves  is  known  prior  to  the  arrival  of  the  vessel.  The  locations  where 
discharge  containers  are  going  to  be  stored  are  not  given.  So,  we  randomly 
assign  areas  of  preference  for  discharge  containers. 


Internal  Trucks 


□ 

— 

Vessel 

— 

Quay  Cranes 


Figure  9.  Yard  Layout 


We  take  60  containers  as  a  base  case  for  our  study.  30  of  those  60 
containers  are  discharge  and  30  are  load  containers.  First,  we  run  the  model  for 
60  containers.  The  ILP  contains  about  1,400,000  equations,  507,000  continuous 
and  487,000  binary  variables.  After  about  28  hours,  the  model  runs  out  of 
memory.  The  lower  bound  it  provides  is  about  336  and  the  best  integer  value  is 
520.  This  means  a  35.38%  relative  gap  (ratio  of  the  difference  between  the 
integer  value  and  the  best  possible).  This  result  and  the  length  of  the 
computation  time  lead  us  to  use  an  optimization  based  heuristic  method  to  get  a 
better  solution  in  less  than  one  hour.  We  use  a  Cascade  Method  (Baker  and 
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Rosenthal  [1998]  and  Baker  et  al.  [2002]).  With  this  method,  we  solve  the  sixty- 
container  problem  part  by  part.  Namely,  we  run  the  model  for  a  small  part  of  data 
at  a  relative  gap  of  5%  for  each  part  and  then  we  run  the  model  again  for  a 
relatively  larger  number  of  containers  after  fixing  some  of  the  variables  to  the 
values  that  we  obtain  at  the  first  run.  This  process  goes  on  until  a  solution  is 
obtained  for  60  containers.  Note  that  each  part  of  the  problem  is  solved  by  an 
exact  algorithm  (CPLEX  Version  10.0.1). 

We  set  up  different  cascade  scenarios  and  run  the  model  for  each  of 
them.  Table  1  shows  the  scenarios  that  we  tested. 


Scenario 

1st 

Step 

2nd 

Step 

3rd 

Step 

4th 

Step 

5th 

Step 

MIP 

Solution 

Run  Time 
(min.) 

1st  Scenario 

Fixed 

— 

0-10 

0-20 

0-30 

0-40 

526 

30 

Optimized 

20 

30 

40 

50 

60 

2nd  Scenario 

Fixed 

— 

0-10 

0-20 

0-30 

479 

41 

Optimized 

20 

30 

50 

60 

3rd  Scenario 

Fixed 

— 

0-10 

0-20 

0-30 

0-40 

486 

52 

Optimized 

25 

35 

45 

55 

60 

4th  Scenario 

Fixed 

— 

0-20 

0-40 

479 

37 

Optimized 

25 

45 

60 

5th  Scenario 

Fixed 

— 

0-25 

528 

86 

Optimized 

25 

60 

Table  1 .  MIP  Solutions  and  Computation  Times  with  Thirty-six  Areas 

In  the  first  scenario,  at  the  first  step,  we  run  the  model  for  20  containers. 
At  the  second  step,  we  fix  the  variables  associated  with  the  first  10  containers  to 
the  values  that  we  obtain  at  the  first  step  and  run  the  model  for  40  containers.  At 
the  third  step,  we  fix  the  variables  associated  with  the  first  20  containers  to  the 
values  that  we  obtain  at  the  second  step.  We  follow  this  process  until  we  reach  a 
solution  for  60  containers.  At  the  fifth  step,  we  run  the  model  for  60  containers 

after  fixing  the  variables  associated  with  the  first  40  variables.  We  get  an  integer 
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solution  with  the  value  of  526  with  this  scenario  in  a  total  of  30  minutes.  The 
fourth  scenario  provides  the  best  result  regarding  the  integer  solution  and  the 
computation  time.  It  provides  an  integer  solution  with  the  value  of  479  in  37 
minutes.  The  relative  gap  is  29.85%  (  better  than  the  35.38%  that  we  got  before). 
The  computation  time  is  also  much  shorter. 

The  assignments  of  RTGs  also  change  with  each  scenario.  Figure  10 
shows  the  assignments  of  RTGs  to  the  areas. 
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Figure  10.  RTG  Assignments  With  36  Areas 


In  Figure  10,  rectangles  represent  areas.  There  are  two  numbers  in  each 

rectangle.  The  number  on  the  right  is  the  area  number.  The  circled  ones  on  the 
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left  show  the  RTG  that  is  assigned  to  that  area.  Shaded  rectangles  represent  the 
areas  where  RTGs  are  initially  located.  Some  of  the  areas  are  not  covered 
because  there  are  several  areas  where  a  discharge  container  can  be  stored. 
Also,  a  load  container  can  usually  be  picked  up  from  multiple  areas. 

We  also  set  up  the  similar  scenarios  for  a  lesser  number  of  areas.  We 
decrease  the  number  of  areas  from  36  to  18,  and  define  new  locations  for  load 
containers  and  initial  locations  for  RTGs.  We  keep  the  rest  of  the  data  the  same 
and  run  the  model  again.  First,  we  run  the  model  for  the  base  case  scenario, 
which  has  no  fixed  values.  The  ILP  contains  about  360,000  equations,  137,000 
continuous  and  127,000  binary  variables.  After  about  28  hours,  the  model  runs 
out  of  memory.  The  lower  bound  it  provides  is  225  and  the  best  integer  value  is 
283.  This  equals  a  20.15%  relative  gap.  Table  2  shows  the  results  of  different 
scenarios  when  we  use  the  cascade  method. 


Scenario 

1st 

Step 

2nd 

Step 

3rd 

Step 

4th 

Step 

5th 

Step 

MIP 

Solution 

Run  Time 
(min.) 

1st  Scenario 

Fixed 

— 

0-10 

0-20 

0-30 

0-40 

264 

3 

Optimized 

20 

30 

40 

50 

60 

2nd  Scenario 

Fixed 

— 

0-10 

0-20 

0-30 

321 

3 

Optimized 

20 

30 

50 

60 

3rd  Scenario 

Fixed 

— 

0-10 

0-20 

0-30 

0-40 

267 

4 

Optimized 

25 

35 

45 

55 

60 

4th  Scenario 

Fixed 

— 

0-20 

0-40 

303 

3.5 

Optimized 

25 

45 

60 

5th  Scenario 

Fixed 

— 

0-25 

344 

5 

Optimized 

25 

60 

Table  2.  MIP  Solutions  and  Computation  Times  with  Eighteen  Areas 


This  time  the  first  scenario  provides  the  best  integer  value.  The  relative 
gap  associated  with  that  integer  value  is  14.77%.  The  computation  time  is  only  3 
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minutes.  This  is  a  reasonable  time  when  compared  to  the  computation  time  of 
the  base  case,  which  is  more  than  28  hours. 

Figure  11  shows  the  assignments  of  RTGs  to  the  areas,  for  the  scenarios 
of  18  areas. 
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Figure  11.  RTG  Assignments  With  18  Areas 

Figure  11  shows  that  the  model  tends  to  assign  the  RTGs  to  adjacent 
areas.  This  is  due  to  the  low  transfer  cost  between  those  areas.  However,  in 
some  situations,  an  RTG  may  be  assigned  to  an  area  that  is  far  away  from  its 
current  location  (e.g.,  in  the  2nd  scenario  RTG5  is  assigned  to  area  01A).  This 
happens  when  there  is  a  container  to  be  handled  in  that  area  and  the  closest 
RTG  is  not  available.  Then,  among  the  available  RTGs,  the  RTG  that  is  located 
in  an  area  associated  with  the  lowest  transfer  cost  is  assigned  to  that  area. 
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V.  CONCLUSION 


This  thesis  is  the  first  to  combine  the  container  handling  and  yard  crane 
deployment  strategies  to  improve  the  efficiency  of  yard  management  in  container 
terminals. 

We  establish  a  model  that  finds  the  best  area  for  a  container  that  is 
unloaded  from  a  vessel,  selects  a  container  to  load  among  several  possible 
areas  and  assigns  yard  cranes  to  the  areas  where  containers  are  stored.  We 
formulate  the  problem  as  an  ILP  and  test  it  with  real  world  data.  Due  to  the 
complexity  of  the  problem,  we  break  the  data  into  small  parts.  Even  for  the 
moderate-size  data,  it  takes  more  than  twenty-four  hours  for  the  ILP  to  provide  a 
near  optimal  solution.  In  order  to  decrease  the  computation  time  and  to  satisfy 
memory  requirement,  we  apply  a  cascade  method  that  solves  the  problem  as  a 
sequence  of  restricted  subproblems.  Each  subproblem  is  restricted  to  a 
sequence  of  containers  and  the  output  of  each  subproblem  provides  an  input  to 
the  next  subproblem.  This  method  provides  better  solutions  than  the  solution  that 
we  get  by  solving  the  problem  directly.  The  cascade  method  also  decreases  the 
computation  time  significantly.  The  results  demonstrate  the  ability  to  combine 
container  handling  and  yard  crane  deployment  in  a  single  model  and  they  verify 
that  cascade  method  works  well  with  the  ILP.  As  a  general  conclusion,  this  model 
might  be  used  to  provide  a  benchmark  for  the  quality  of  the  solutions  provided  by 
heuristic  algorithms. 

In  this  study,  we  combine  two  of  the  three  major  operations  in  a  container 
terminal.  Further  studies  may  add  the  landside  operations  and  provide  more 
general  solutions  for  container  terminal  operations.  Additionally,  heuristic 
algorithms  may  be  used  to  provide  real-time  solutions. 
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